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ABSTRACT 



Aims. A model to describe cosmic ray spectra in the energy region from 10 to 10 eV is suggested based on the assumption 
that Galactic cosmic ray flux is a mixture of fluxes accelerated by shocks from nova and supernova of different types. 
Methods. We analyze recent experimental data on cosmic ray spectra obtained in direct measurements above the atmosphere 
and data obtained with ground Extensive Air Shower arrays. 

Results. The model of the three classes of cosmic ray sources is consistent with direct experimental data on cosmic ray 
elemental spectra and gives a smooth transition from the all particle spectrum measured in the direct experiments to the all 
particle spectrum measured with EAS. 



Key words, cosmic rays - ISM:supernova remnants - ISM:abundances 



1. Introduction 

It is thought that cosmic rays below 10 17 eV are acceler- 
ated by shocks of supernova explosions, with all compo- 
nents having the same rigidity spectra. However, recent 
experimental data show that spectral indices of elemen- 
tal spectra are different. Therefore we suggest that the 
concepti of a single population of sources with the same 
index of rigidity spectra for all elements is inaccurate We 
describe the experimental data using the idea that the 
measured cosmic ray flux is generated by three different 
classes of sources. It is suggested that each class of sources 
generates the spectrum that is power-law by rigidity with 
its specific spectral index and maximal rigidity. We ana- 
lyze cosmic ray spectra, measured by direct methods, and 
all-particle spectra, measured in the extensive air showers 
(EAS). The data used are shown in Table 1. 



2. Basis for the model 

The model is based on the preliminary experimental data 
obtained in the ATIC-2 experiment (Wefel et al, I2UU51 
Panov et al, I2HB51 Batkov et al, |2TM|) that fill the gap 
between measurements with magnetic spectrometers be- 
low 100 GeV (Haino et al, I2UU1 Alcaraz et al, I2UUU1 
Boezio et al, 120030 and emulsion experiments above 10 
TeV (Zatsepin et al,EM| Asakimori et al. 119981 Derbina 
et al. I2005f) . The proton and helium spectra obtained in 
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Fig. 1. Proton and He spectra. Dashed lines are described 
in Section 3, solid lines are described in Section 5. 



these direct measurements are shown in Fig.l. The prelim- 
inary ATIC-2 data show that the proton spectrum obeys 
a power law with a spectral index of 7 = 2.63 ± 0.01 in 
the energy region between 3 x 10 2 and 10 4 GeV . However, 
the extrapolation of this spectrum to higher energy does 
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Table 1. Experiments used in this paper 



Experiment 


Technique 


Site 


Reference 


AMS 


Magnetic spectrometer 


Spacecraft 


Alcaraz et al., 2000 


CAPRICE 


Magnetic spectrometer 


Balloon 


Boezio et aL.ISUUSl 


BESS-TEV 


Magnetic spectrometer 


Balloon 


Haino et al.,|2T5u3| 


ATIC-2 


Calorimeter 


Balloon 


Wefel et al.. 120051 Panov et al.. 120031 


SOKOL 


Calorimeter 


Spacecraft 


Ivanenko et al.. 119931 


JACEE 


Emulsion chamber 


Balloon 


Asakimon et al..|1998| Takahashi et al..|1998| 


1\/TT TRT?T? 
1VL U rJrjrj 


Emulsion chamber 


Balloon 


ZjatseDin et ai.iiyyoi iiyy4i 


RUNJOB 


Emulsion chamber 


Balloon 


Derbina et al.,|2T503| 


HEAO 


Cerenkov counter 


Spacecraft 


Eneelmann et al. .119901 


CRN 


Transition radiation 


Spacecraft 


Miiller et al.,[TMT| 


TRACER 


Transition radiation 


Balloon 


Miiller et al.,|5003| 


TIC 


Calorimeter 


Balloon 


Adams et al.,|1997| 


KASCADE 


EAS 


Ground based 


Roth et al.,|2003| 


HEGRA-AIROBIC 


EAS+Cerenkov light 


Ground based 


Arqueros et al.|2000| 


CASA-BLANCA 


EAS+ Cerenkov light 


Ground based 


Fowler al.,|2000| 


DICE 


Cerenkov light 


Ground based 


Kieda et al-ITOOOl 


TUNKA 


Cerenkov light 


Ground based 


Budnev et al.. 120051 



not agree with the emulsion experiment data. Let us sup- 
pose that this difference is real, and that we see the end of 
acceleration for the first class of sources, while sources of 
the second class accelerate cosmic rays efficiently up to the 
knee region at about 3 x 10 15 eV. Initially we suppose that 
the steeper spectrum below 3 x 10 2 GeV is connected with 
weak reacceleration of cosmic rays during their propaga- 
tion from the sources to the Earth. The helium spectrum 
is flatter than the proton one. This difference may be ex- 
plained by the assumption that the second type of sources 
has a different spectral index, and that the elemental com- 
position for these two classes of sources arc different. 

3. Description of the model 

The equilibrium spectra in the Galaxy will be described 
in the framework of the Leaky Box model: 



I(R) = 



Q(R) x Tesc (R) 
1 + A esc (i?)/ Xp 



(1) 



where R is the particle rigidity, Q(R) is the source spec- 
trum, T esc is the cosmic ray life time for the escape, A p 
is the particle mean free path for interaction, and A esc is 
the escape length (A esc (i?) = p x v x T esc (R), where p is 
the density of the interstellar medium and v is the particle 
velocity). 

The Leaky Box model without reacceleration ade- 
quately describes the spectra of nuclei and the ratios of 
secondary nuclei to primary ones in the HEAO-3-C2 ex- 
periment below 35 GeV/nucleon (Engelmann et al'1990), 
if 

A esc (i?) = 34.1 x R-° 6 g/cm 2 , for R > = 4.4 GV (2) 

\ esc {R) = 10.8 x /?, g/cm 2 for R < 4.4 GV 
and Q(R) ~ iT 2 23 ; 



but we assume that at high energy X esc (R) ~ Rr 1 ' 3 , that 
corresponds to the Kolmogorov type of magnetic turbu- 
lence, and that the dependence in eq.(2) at R < 100 GV 
results from reacceleration of cosmic rays during their 
propagation in the Galaxy. According to Osborne and 
Ptuskin ( 198SJ) reacceleration may be taken into account 
as follows: 

X esc (R) = 4.2 x (R/Ro)- 1 / 3 x [1 + (R/R Q )- 2 / 3 } g/cm 2 , 

(3) 

where i? = 5.5 GV. It is assumed that 



Q(R) - R- a x cj>(R), 



(4) 



where a is the index of the source spectrum, and the func- 
tion 4>{R) describes smooth transition from the spectral 
index in the region of effective acceleration for each type 
of source to the spectral index after termination of this 
process. 

4>(R) = ^ + (R/R max ) 2 ] (l ~ lk)l2 (5) 

where 7 ~ a + 0.33 is the spectral index in the region of ef- 
fective acceleration (at high enough energy), and 7^ is the 
spectral index after termination of effective acceleration. 
It is assumed that spectra are simple power laws with the 
index 7^ after termination of effective acceleration. 

As later we will fit spectra of nuclear groups and the all 
particle spectrum, we convert rigidity spectra to spectra 
by energy per particle E. 



1(E) 



Q P (R) x r esc (R) dR 
l + \ esc {R)/\ p X dE 



R = — x ^JE 2 + 2m p x Ax E 



dR 



E + m p A 



(G) 



(7) 



dE Z y/E 2 + 2m p x Ax E 

where Z and A are particle charge and atomic weight, and 
m p is proton mass. 
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Fig. 2. Spectra of nuclear groups. Dashed lines are de- 
scribed in Section 3, solid lines are described in Section 
5. 



The spectrum of each cosmic ray nuclear group is the 
sum of spectra from the different classes of sources. 

We described sources of class I with values of a — 2.3 
and Kmax — 50 TV, and sources of class II with a = 2.1 
and Rmax — 4 PV. The intensities for various cosmic ray 
groups were chosen to fit both the data of direct measure- 
ments and the data of EAS in the high energy region. The 
predictions of the model along with experimental data for 
cosmic ray groups and the all particle spectrum are shown 
in Figures 1-3 with dashed lines (above 10 4 GeV this line 
follows the solid line, which is described in Section 5). The 
parameters of the model are shown in Tables 1 and 2. 

Above 300 GeV the proton spectrum (Fig.l) is fitted 
closely by the model. The steepening of the helium spec- 
trum above 10 TeV is not as clear as for proton one, but 
the line does not contradict the experimental data. The 
RUNJOB data for helium are much lower than the data of 
ATIC-2 and JACEE. This also holds true for other groups 
of nuclei. In the region below ~ 300 GeV the model does 
not agree with the experimental data. An identical be- 
havior is seen for the spectra of heavier nuclear groups 
shown in Fig. 2 and for the all-particle spectrum shown 
in Fig. 3. We will discuss this discrepancy in Section 5. 
The parameters for class II that determine the knee re- 
gion were chosen to fit the experimental data of HEGRA- 
AIROBICC, KASCADE and TUNKA, while the data of 
CASA-BLANKA and DICE appear to be below the model 
line. 



ATIC-2 
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Fig. 3. All particle spectrum. The dashed line is described 
in Section 3, the solid line is described in Section 5. 



4. Interpretation of the model 

Our analysis is not the first attempt to link data of direct 
measurements with data of EAS and to fit the knee. Some 
previous models arc phcnomenological, e.g. the poligonato 
model of Horandel ( 2003 ) and some are based on physical 
grounds. In a number of papers, the origin of the knee 
is explained by the structure of magnetic fields in the 
Galaxy, or the change of inelastic interaction character- 
istics above several PeV. Erlykin and Wolfendale l|1997[) 
assumed that the knee arises from a single nearby cosmic 
ray source. A detailed review of various models was made 

mm. 



by Horandel 

As two classes of cosmic ray sources are needed for 
our model, we use the model of Bierman ( I1993f) where 
two types of astrophysical objects are supposed as sources 
of the Galactic cosmic rays. The first type is stars with 
masses from 8 to 15 of solar mass that exploded into the 
interstellar medium (ISM). The expected chemical com- 
position generated by these objects is close to the compo- 
sition of the ISM. The second type is massive stars explod- 
ing into their own stellar wind. The stars of 15 - 25 solar 
masses known as Red Supergiants (RSG) possess a mod- 
erate stellar wind and the stars of mass greater than 25 so- 
lar masses known as Wolf-Rayet (WR) stars have a strong 
stellar wind. The wind stars are the late stage of evolu- 
tion of massive O and B stars, which arc clustered into 
OB associations. The combined effect of collective stellar 
winds and supernovae within OB associations leads to the 
production of a superbubble (Heiles 119791 Bruhweiler et 
al,rT980, Tomisaka ll992fl . i.e. a huge (several hundred par- 
sec) region of hot (~ 10 6 K), low-density (~ 10~ 3 cm~ 3 ) 
gas, surrounded by a cold neutral hydrogen supershell. 
It was supposed that the Sun is located within such a su- 
perbubble (Kafatos et al.. H981f) . It was shown (Kafatos et 
al. rtWn Streitmatter et al.. 1X5551 Lingenfelter et al. 1217001 
Higdon and Lingenfelter, [2003 Higdon and Lingenfelter, 
2005, Streitmatter and JonespOOS) that in this case many 
observed features of cosmic rays can be explained (age, 
anisotropy, some features of isotopic composition) . We as- 
sume that our first class comes from with the explosions 
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of isolated stars into the ISM, and that our second class 
comes from supernovae within the local superbubble. 

We assume that the second type of source generates 
cosmic rays with natter spectra and with higher maximal 
energy. These are those cosmic rays that are studied by 
ground EAS arrays, and that produce the knee in the PeV 
energy region. 

5. Energy spectra below ~ 300 GeV per nucleon 

We now direct our attention to the energy region below 
~ 300 GeV per nucleon. As one can see from Fig.l - 3 
(dashed lines), the model with two classes of sources does 
not fit the preliminary ATIC-2 experimental data well. 
The weak reacceleration is unable to fit the spectra in the 
low energy region and it may be suggested that most of 
cosmic rays below 300 GeV per nucleon are accelerated 
in sources of a third class. The parameters for sources 
of the third class required to fit the experimental data 
are the following: a = 2.57; R max = 200 GV; j k = 4.5. 
In the region of such low rigidities, the solar modulation 
must be taken into account. We have used the eq. from 
Boczio( 1997) for this purpose with modulation parameter 
$ = 600 MV. 



flux mod {E) = flux(E + Ze x $) x P, 
^_ E 2 + 2m p E 



(8) 



{E + Zex $) 2 + 2m p (E + Ze x $) ' 

where e is the charge of the electron, m p is the proton 
mass and E is the kinetic energy of the nucleus. 

The intensities for various cosmic ray components were 
chosen to fit the experimental data in Figures 1-3. The 
fit in the version of model with three classes of sources is 
shown in Fig. 1-3 with solid lines. 

To check our model in this low energy region we 
plot experimental data for separate heavy elements, oxy- 
gen and iron, obtained from HEAO-3-C2 (Engelmann et 
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Fig. 5. Top panel: Fluxes of cosmic ray groups in the 
model. Line numbers 1: protons, 2: Helium, 3: CNO group, 
4: Ne - S group, 5: Fe group (Z> 17). Bottom panel: 
< InA > vs energy; small squares: preliminary ATIC-2, 
large solid squares: JACEE, open squares: RUNJOB 



al., HMTjl . CRN (Muller et al., ITMTjl . and preliminary 
TRACER results (Muller et al., |2005) and the prelimi- 
nary ATIC-2 data (Panov et al.. I2()()5|l along with two 
versions of the model in Fig. 4. The model of three classes 
of sources fits the experimental data well. The model also 
fits other individual spectra. 

Nova stars are candidates for the third class of sources. 
During nova explosion, an expanding shell is produced, 
and the broad features of nova explosion are similar to 
these of the explosion of supernovae. Power requirements 
on the power are also satisfied: the mean energy of explo- 
sion is W ~ 10 46 — 10 47 erg., the frequency of nova ex- 
plosions is ~ 100 year -1 ( Pskovskv. l2()()()|l . So, the power 
available for cosmic ray generation is: 



(10 4b - 10 47 ) x 100 
3 x 10 7 



3 x 10 40 - 3 x 10 41 erg/sec 



It is not lower than the power supplied by supernova. 
Thus, this source class could generate cosmic rays with 
a steeper spectrum (a = 2.57) and of significantly lower 
maximal rigidity ( 200 GV). 
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Table 2. Parameters for three classes of sources. Class 
I: SN into ISM; class II: SN in the Superbubble, class III: 
Novae. 



Class 


a 


Umax [GV] 


7 


7fc 


I 


2.3 


5 x 10 4 


2.63 


8 


II 


2.1 


4 x 10 6 


2.43 


4.5 


III 


2.57 


2 x 10 2 


2.9 


4.5 



Table 3. Flux xE 2 - 75 [m^sec^ster- 1 GeV 175 for cos- 
mic ray groups at E = 10 TeV 





Z eff 


Fluxx£^ Y5 




(A P ) 


class I 


class II 


class III 


H 


1(74) 


1.4 x 10 4 


6.25 x 10 3 


0.6 


He 


2(18) 


8.5 x 10 s 


8.5 x 10 3 


1.5 


CNO 


7(5.8) 


6.75 x 10 3 


1.8 x 10 3 


30 


Ne-S 


12(3.5) 


5.5 x 10 3 


1.5 x 10 3 


110 


Z > 17 


20(2.4) 


3.5 x 10 3 


1.2 x 10 3 


750 



6. Parameters of the model, fluxes of different 
cosmic ray groups and <lnA>. 

The parameters of the model used to fit spectra of cosmic 
ray groups are shown in Table 1 and Table 2. Fig. 5 shows 
fluxes of cosmic ray groups and the mean value of <lnA> 
in the model. The special feature of our model is that 
charge composition is not a monotonic function of energy: 
it becomes 'lighter' before the knee, and only above the 
knee becomes heavier. 



7. Conclusion 

We proceed from the assumption that the difference be- 
tween proton and helium spectral indices found in the 
JACEE experiment above 10 TeV (Asakimori et al !1998|) . 
and supported by the preliminary ATIC-2 analysis (Wefel 
et al.. 12005(1 in the energy region above 100 GeV is real. 
We suppose that this difference results from the fact that 
cosmic ray fluxes are a mixture of fluxes from sources with 
different spectral indices and different maximal energy. We 
checked a simplest model assuming two classes of sources. 
We showed that elemental cosmic ray spectra pointed to 
a 'kink' at the rigidity of about 50 TV, and connected 
this with the assumption that one class of sources termi- 
nates its effective acceleration at this rigidity. This source 
class may be identified with the supernova explosions into 
the ISM. The second source class, presumably supernovae 
within the local superbubble, accelerates cosmic rays up to 
rigidity of 4 PV, producing the classic knee. However, the 
assumption of these two classes of sources is inadequate to 
fit the energy spectra below 300 GeV per nucleon along 
with the spectra in the high energy region, even taking 
into account minimal reacceleration of cosmic rays in the 
Galaxy. We assume that the contribution of nova stars is 
essential in the energy region below ~ 300 GeV per nu- 
cleon. 
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